Pseudogap Induced Antiferromagnetic Spin Correlation 
in High- Temperature Superconductors 
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The pseudogap phenomena observed on cuprate high temperature superconductors are investigated 
based on the exact diagonalization method on the finite cluster t-J model. The results show the 
presence of the gap-like behavior in the temperature dependence of various magnetic properties; 
the NMR relaxation rate, the neutron scattering intensity and the static susceptibility. The cal- 
culated spin correlation function indicates that the pseudogap behavior arises associated with the 
development of the antiferromagnetic spin correlation with decreasing the temperature. The nu- 
merical results are presented to clarify the model parameter dependence, that covers the realistic 
experimental situation. The effect of the next-nearest neighbor hopping t' is also studied. 



I. INTRODUCTION 

The high-temperature superconductors exhibit various 
anomalous metallic properties in their normal phases. A 
gap-like behavior, first discovered by the NMR relaxation 
rate measurement |lj as one of these properties has at- 
tracted a lot of current interest. It is characteristic to the 
underdoped cuprates and observed as a broad peak in the 
X/T\T-T curve at a slightly higher temperature than the 
superconducting transition point T c . The possible rela- 
tion of the pseudogap phenomena with the origin of the 
superconductivity underlies the reason of such interest. 
The behavior has now been detected in the temperature 
dependence of various other physical quantities; the neu- 
tron scattering intensity the magnetic susceptibil- 
ity Q, the resistivity g, the Hall coefficient ||] and the 
angle-resolved photo-emission spectrum (ARPES) |Q. 

A lot of theories have been proposed to explain the 
phenomena. Since they were observed only for com- 
pounds with double Cu02 layers like YB^C^Og for 
the first time, the interlayer coupling was supposed to 
play the predominant role on the gap formation. || p"c| ] 
However, after the discovery of the phenomena on the 
mono-layer oxide HgBa2Cu04 by the NMR measurement 
[ lllfl , most theories are based on the model with a single 
CuC>2 plane, such as the square-lattice t-J model which 
has been extensively studied theoretically in the descrip- 
tion of the cuprate superconductivity. |l2[ Based on the 
concept of the spin-charge separation, | |13| | Tanamoto et 
al. Q showed that the t-J model exhibits a gap- like be- 
havior in its spin excitation spectrum as the result of the 
spinon condensation by applying the mean field approx- 
imation to the resonating valence bond basis. Onoda et 
al. [jl5| also derived the occurrence of the gap-like be- 
havior in the temperature dependence of the resistivity 
within the same model using the slave-boson mean field 
approximation. These results, however, crucially depend 
on the nature of the approximation, the validity of which 
is not so well established. 

Miyake and Narikiyo, on the other hand, suggested 
that the pseudogap rather originates from the enhance- 



ment of the antiferromagnetic spin fluctuation paying at- 
tention on the nesting effect based on the phcnomenologi- 
cal itinerant-localized duality theory. |l6f| In our previous 
study, we also proposed a simple explanation based on 
the well-known property of the Heisenberg magnets. (Tif l 
We were successful in demonstrating the presence of the 
gap-like behavior in the temperature dependence of the 
NMR relaxation rate by numerically diagonalizing the fi- 
nite cluster t-J model for sufficiently large J. Our result 
shows it is interpreted as a crossover from the paramag- 
netic to the magnetically well correlated state induced 
by the rapid development of the antiferromagnetic spin 
correlation. 

Lately lots of interest have been paid on the role of the 
superconducting fluctuations. For instance, Koikegami 
and Yamada |18| , based on the d-p model, indicated that 
superconducting fluctuations will survive at higher tem- 
perature than the transition temperature T c based on the 
fluctuation exchange approximation. This means that 
the gap-like behavior is regarded as a precursor of the 
gap formation of the strong coupling superconductivity 
as was predicted by the self-consistent i-matrix approxi- 
mation by Yanase and Yamada . With the use of the 
t-matrix approximation Kobayashi et al ppf ] showed that 
the gap-like behavior is induced by the appearance of the 
superconducting fluctuations while the antiferromagnetic 
spin fluctuation is suppressed on the contrary. Onoda 
and Imada claim that the competition between the 
superconducting fluctuations and antiferromagnetic spin 
fluctuations is necessary for deriving the gap behavior. 

In spite of lots of theoretical proposals, the final an- 
swer is still far from being established. Most theories, 
however, seem to have an agreement on the important 
role of aniferromagnetic spin correlation on the pseudo- 
gap formation. In the present paper, in order to confirm 
the validity of our mechanism we extend our numerical 
investigation and derive the temperature dependence of 
various physical properties, including \jT\T, in a wider- 
parameter region that covers realistic situations for high- 
T c cuprates. The effect of the next-nearest neighbor hole 
hopping with a coupling t' is also studied. The effect of 
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the antiferromagnetic spin correlation on the pseudogap 
formation is then clarified without resorting to any ap- 
proximation. If we take into account the lack of reliable 
approximation scheme in treating the strongly correlated 
systems, the numerical analysis of the model has its own 
significance in order to check the validity of various ex- 
planations. 

II. MECHANISM OF PSEUDOGAP 



imaginary part of the dynamical spin susceptibility of 
conduction electrons lmx(q,to) and the static spin cor- 
relation function (S?S?). In actual numerical estima- 
tions of Irnx(q, the (5-function is approximated by the 
Lorentzian distribution with some small width. Since we 
need all the eigenvalues and eigenvectors, the cluster size 
of the model was limited by the available disk space of 
the computer. We show in the following the results for 
the V^O x \/T0 cluster under the periodic boundary con- 
dition. 



According to our proposal, the local spin excitations 
are suppressed around a characteristic temperature, of 
the the same magnitude of the antiferromagnetic cou- 
pling J, associated with the growth of the short-range 
antiferromagnetic correlation as we decrease the tem- 
perature. As a result various magnetic quantities will 
show their particular T-dependence in reflecting the re- 
duced local excitations. In fact the high temperature 
series expansion for the square-lattice antiferromagnetic 
Heisenberg model shows a broad peak in the temperature 
dependence of the NMR relaxation rate \fT\T around 
T ~ J. [0 Since the hole motion acts to destroy the 
short-range antiferromagnetic order, the hole doping will 
reduce the crossover temperature in agreement with the 
observed hole doping dependence of the pseudogap tem- 
perature. Thus the pseudogap behavior is simply re- 
alized as the manifestation of the growth of the short- 
range antiferromagnetic correlation. The above picture 
has been supported by our numerical exact diagonaliza- 
tion study on the two-dimensional t-J model with 10 lat- 
tice sites. [|l7| We found that for a small hole concentra- 
tion (6 = 0.1) the T dependence of X/T\T does show a 
broad peak around the temperature where the nearest 
neighbor antiferromagnetic spin correlation shows rapid 
growing, for a sufficiently large antiferromagnetic cou- 
pling (J/t = 0.6). In the next section we perform the 
same calculation to investigate how our previous results 
are modified by using realistic parameters for the cuprate 
superconductors and also to examine the presence of pos- 
sible gap-like behaviors for various magnetic properties. 



III. NUMERICAL ANALYSIS OF T-J MODEL 

Let us first consider the standard square-lattice t-J 
model Hamiltonian defined by 

H = -*E<ij>, < r(c j > Ci ^ + 4,a c i,cr) 

+ ^E< £ j>(Si-S j -in i n j ), (1) 

where t is the nearest neighbor electron hopping integral 
and J is the antiferromagnetic Heisenberg exchange con- 
stant between spins on adjacent lattice sites. Through- 
out the paper, all the energies are measured in units of 
t. With the use of the calculated eigenvalues and eigen- 
vectors, we evaluated the temperature dependence of the 



A. NMR relaxation rate 

In order to obtain the temperature dependence of the 
NMR relaxation rate, we evaluate \jT\T by the following 
formula: 



— — = lira - Vlmyffl.ij). 



(2) 



The effect of the q-dependence of the hyperfine form fac- 
tor is neglected, for simplicity. To clarify the role of an- 
tiferromagnetic spin correlation, we have also evaluated 
the T dependence of the Q = (n, n) component of the 
spin correlation function, 



(j-i)-(*+y) 



(S?S? 



(3) 



In our previous study |L7| the following nearest-neighbor 
spin correlation function was used instead for the pur- 
pose, 



^ = ^Ei E 



+ P i- 



(4) 



p=±x,±y 



We are using S{Q) here to emphasize the continuity that 
there exists the antiferromagnetic order in the limit of 
no doping. Note that the terms, magnetically singlet and 
antiferromagnetic, are used in nearly the same meaning. 
We are only concerned with the short range magnetic 
correlation of the system of finite clusters. According 
to the above definitions there is no reason to make defi- 
nite distinction between S(Q) and C\. Actually both of 
S(Q) and C\ show almost no essential difference in their 
temperature dependence p3| . 

We show in Fig. results of T dependence of (a) \JT X T 
and (b) S(Q) for the undoped case (S = 0) for J=0.6, 0.5, 
and 0.4. Only for J=0.4, the temperature dependence of 
S(Q) has already been reported by Tohyama et al. p4| ] 
For cuprate superconductors the value J = 0.4 is esti- 
mated experimentally. All the calculated T dependence 
of \jT\T exhibits a broad peak around T ~ J where 
S(Q) shows the significant increase. The above behav- 
ior of \jT\T for S = is, for instance, consistent with 
the high temperature series expansion for the Heisen- 
berg model. |£2 25 The quantum Monte Carlo simulation 
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study by Sandvik and Scalapino |p6[ , however, showed 
no peak behavior in the \jT\T-T curve. No peak is 
also observed by the NMR measurement for the undoped 
La2CuC>4. [^7j The reason is the temperature in these 
studies is not so high enough to cover the crossover re- 
gion of the order of J. 

The same results are shown in Figs.|](a) and ||(b) for 
the one-hole case (<5 = 0.1). Peaks in the T dependence 
of 1/TiT for J=0.5 and 0.6 are shifted to lower tem- 
perature where the rapid growth of S(Q) is observed, in 
agreement with our mechanism. Although no peak ap- 
pears for the case J=0.4, a slight hump still exists around 
the same temperature. Such a behavior will be detected 
as a deviation from the Curie- Weiss like temperature de- 
pendence of 1/T\T, rather than a peak as was actually 
observed for La2- x Sr a; Cu04. p8| The same calculation 



for the hole-doped 4x4 t- J cluster 29 does not show any 
gap-like anomalies. This is probably because the temper- 
ature range didn't cover the region where the pseudogap 
is expected to appear. The present results indicate that 
the gap-like behavior appears at much higher tempera- 
ture than the calculated lowest excitation gap A of the 
cluster. For example, the broad peak appears around 
T ~ 0.4, while the estimated A is 0.146, for J = 0.6 and 
(5=0.1. |l7]] Therefore it excludes the possibility that the 
pseudogap is caused by the discreteness of energy levels 
of finite systems. The same results for the two-hole sys- 
tem (S — 0.2) are shown in Fig s|(a) and |(b). Although 
no gap-like behavior is observed, they are also consistent 
with our view because of the absence of any significant 
enhancement of the antiferromagnetic spin correlation. 
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FIG. 1. Temperature dependence of (a) 1/TiT and (b) 
S(Q) for 5 = 0. 
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FIG. 2. Temperature dependence of (a) \jT\T and (b) 
S(Q) for (5 = 0.1. 
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FIG. 3. Temperature dependence of (a) 1/TiT and (b) 
S(Q) for (5 = 0.2. 



B. Neutron scattering intensity 

The measurement of the neutron scattering intensity 
of YBa2Cu306.6 by Sternlieb et al. ||| indicated the pres- 
ence of the pseudogap behavior in the temperature de- 
pendence of the g-integrated dynamical susceptibility 



dqlmx(q,uj)- 



(5) 



We show, in Figs^j, || and [| the temperature depen- 
dence of (a) Imx(Q,uj) and (b) Imy(w) for a small to 
(=0.01t) for (5=0, 0.1 and 0.2. Figsjfa), |(a) and |(a) 
show the monotonic increase of Imx(Q, w) for all the val- 
ues of J with decreasing T, independent of the doping 
<5. On the other hand, g-integrated local susceptibility 
in Figs.^(b), ||(b) and ||(b) reveals a gap- like behavior 
around the same crossover temperature as that of \jT\T 
for (5 < 0.1. These results are in good agreement with the 
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above neutron scattering experiment. As for the former, 
the low-frequency limit of lmx(Q, u>) is given by 



x(Q,o) 



cx X 2 (g,o), (w-fO), 



where Tq is the Lorentzian distribution width of the fre- 
quency spectrum. The monotonic increase of lmx(Q,t-<j) 
with decreasing the temperature is therefore consistent 
with the absence of the pseudogap behavior in the In- 
dependence of UT 2G , i-e. T 2 2 G oc x _1 (Q,0), by the NMR 
measurement. BOI] 
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FIG. 4. Temperature dependence of (a) ImxiQ,^) and 
(b) ImxO) for (5 = 0. 
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FIG. 5. Temperature dependence of (a) Imx(Q,a>) and 
(b) ImxH for S = 0.1. 
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FIG. 6. Temperature dependence of (a) Imx(Q,w) and 
(b) ImxH for 5 = 0.2. 



C. Magnetic susceptibility 

The pseudogap behavior has also been observed in the 
temperature dependence of the static magnetic suscepti- 
bility x- 01 The calculated results of x f° r <S=0, 0.1 and 
0.2 are shown in Figs. 0, || and [|, respectively. In all 
the figures we can see the crossover behaviors around the 
same temperature as those of \jT\T and Imx(w) in the 
underdoped case (S < 0.1), even for the realistic value 
J = 0.4. As shown in Fig. 9, peaks appear at lower 
temperature for <5 = 0.2. However they result from the 
different origin. If there are two holes in a finite cluster 
model, its ground state always becomes singlet and the 
total energy has some small excitation gap. Therefore 
the susceptibility shows a peak around the temperature 
corresponding to the lowest excitation gap of the system 
that comes from these finite size effects. The behavior 
has nothing to do with our mechanism. (The same cal- 
culation of x f° r J = 0.4 by Tohyama et al. yielded 
a peak even for 6 = 0.3.) 
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FIG. 7. Temperature dependence of x for 8 — 0. 
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FIG. 8. Temperature dependence of y for <5 = 0.1. 
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FIG. 9. Temperature dependence of y for 8 = 0.2. 

To conclude, all the results presented above support 
our view that the pseudogap behavior observed for vari- 
ous magnetic properties are induced by the growth of the 
antifcrromagnetic spin correlation. 

IV. NEXT-NEAREST NEIGHBOR HOPPING T' 

In the previous section we showed that the antiferro- 
magnetic spin correlation produces a gap-like behavior in 
the T-dependence of 1/T\T, Imy(o;) and \ around a com- 
mon crossover temperature. The characteristic tempera- 
ture for S = 0.1 is roughly estimated as ~ J/2, as shown 
in Figs. |^(a), ||(b) and ||. It is, however, much higher 
than observed ones for real cuprates. This is partly be- 
cause the effect of the hole hopping t is overestimated 
compared with that of J due to the finite size effect. As 
another possible reason, let us discuss below the effect of 
including the next-nearest neighbor hole hopping, 



<i,j>',<T 



into the t-J Hamiltonian ([!]), where 2<ij>' means the 
summation over all the next-nearest sites. The term is 
originally introduced to explain the dynamics of the cor- 
related motion of spin and hole-hopping degrees of free- 
dom in some cuprates. The value t' ~ —0.3 is estimated 
to be suitable for Sr2Cu02Cl2 by the exact diagonaliza- 
tion study. To clarify the effect of f ' on the pseudogap 
formation, we show the T dependence of \jT\T for var- 
ious values of t' in Figs. [H](a) for J = 0.6 and [To](b) for 
J = 0.4 (5 = 0.1). From Fig. |o|(a) we can see that the 
peak temperatures are lowered by the effect. The nar- 
rowing of the peak width takes place at the same time. 
As the result the slight anomaly in Fig.^(a) for J = 0.4, 
shifted to lower temperature, now becomes evident if we 
assume t' — —0.3 as shown Fig.[To|(b). Therefore it is 
possible to predict the crossover temperature around a 
realistic temperature range ~ J/4 by including the effect 
of the next-nearest neighbor hopping. 

As we increases the value of t', there appears, however, 
another difficulty. Because the peak temperature is low- 
ered and becomes comparable with the order of the lowest 
excitation gap of the cluster, it is very difficult to distin- 
guish the peak behavior from the one that comes from 
the discreteness of the energy levels of the system. Fur- 
ther calculations on larger size clusters will be necessary 
to clarify the situation and to obtain reliable conclusions. 
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FIG. 10. Temperature dependence of 1/TiT for 5 = 0.1 
in the case of (a) J = 0.6 and (b) J = 0.4. 



V. DISCUSSION 

The results of the present study clearly showed the 
presence of gap-like behaviors in the temperature de- 
pendence of various magnetic properties around a single 
crossover temperature due to the development of anti- 
ferromagnetic spin correlation. Experimentally detailed 
gap-like behaviors differ slightly with each other depend- 
ing on materials. In the case of La2- a; Sr 2 ,Cu04, for in- 
stance, the crossover temperature of y (T*) Q is slightly 
higher than that of 1/T X T (T* MR ) Q. The recent NMR 
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Knight shift measurement for Bi2Sr2CaCu20s+<5 also in- 
dicates the presence of two different crossover tempera- 
tures. 32 Stimulated by these discoveries, it is argued 
that the higher crossover temperature (the large pseudo- 
gap) will be associated with the growth of the antiferro- 
magnetic spin correlation. On the other hand, the lower 
one (the small pseudogap) will be due to the singlet (or 
superconducting pair) formation, i.e. T* will be assigned 
to the former, while ?n MR the latter. From our approach 
based on the finite cluster model, it will be difficult to 
predict the presence of double pseudogap phenomena. 

Our model Hamiltonian employed here is quite simpli- 
fied one. In order to deal with actual systems we have 
to include various additional interactions specific to each 
system. In the present paper we are particularly inter- 
ested in properties that are common to all the cuprate 
superconductors. Therefore in order to compare with ex- 
periments, we have to be aware that the observed proper- 
ties are intrinsic to all of them or not. It will also be im- 
portant to confirm whether the presence of two crossover 
temperatures is common to all the cuprates systems. 

As was pointed out in the Introduction, spin-gap phe- 
nomena is characteristic to the underdoped cuprates 
where the superconducting transition temperature T c de- 
creases proportional to the hole doping concentration. 
On the other hand the pseudogap temperature increases. 
Because of this opposite doping dependence, it seems to 
be difficult to associate the pseudo-gap phenomena with 
the superconducting pairing fluctuations. The antiferro- 
magnetic correlation, therefore, plays the dominant role 
on reducing the low energy excitations of the system as 
far as the magnetic freedoms are concerned. Only around 
the critical temperature T c , they will be further modified 
by the presence of the superconducting transition. These 
effect are limited to the very low energy regions. The 
present mechanism is based on the well-known proper- 
ties of low-dimensional Heisenberg magnets. The overall 
doping dependence of the crossover temperature is also 
in agreements with experiments. It decreases with in- 
creasing hole doping, whereas the transition temperature 
T c increases towards the optimum condition. Because of 
these reasons we suppose that our mechanism is one of 
most probable candidates for the phenomena. 

Because the model treated here is a single-hole system 
(S = 0.1 for the 10 site cluster), its ground state is dou- 
blet. Even in this doublet case we can derive the gap- like 
behavior in the temperature dependence of \jT{T . If the 
same calculation is done on a two-hole system of a larger 
cluster with the same doping 6, a slightly different prop- 
erties will be obtained because its ground state now be- 
comes singlet. We expect that the gap-like behavior will 
then be more evident. The difference becomes less con- 
spicuous, if we increase the size of clusters, for it results 
from the size effect. Though it is very difficult because 
of various limitations concerning the available resources 
of computer systems, it will be interesting to extend our 
present numerical diagonalization study to the two-hole 
system with larger cluster size to clarified the situation. 



VI. SUMMARY 

We have shown that the pseudogap behaviors ob- 
served in the high-temperature superconductors originate 
from the growth of the antiferromagnetic spin correlation 
based on the exact diagonalization study on the finite 
cluster t-J model in a realistic parameter region for the 
real cuprates. We also found that the pseudogap temper- 
ature is lowered by including the next-nearest neighbor 
hopping effect. 
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